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Abstract. We show that the generalized Bell inequality is violated in the 
extended Heisenberg model when the temperature is below a threshold value. 
The threshold temperature values are obtained by constructing exact solutions of 
the model using the temperature-dependent correlation functions. The effect due 
to the presence of external magnetic field is also illustrated. 
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1. Introduction 

An intriguing aspect of quantum mechanics is the lack of a local realistic description 
that could reproduce the necessary correlations for the experimental outcomes in 
composite systems PQ. This lack of local realism can be investigated using the 
entangled state as discussed in the original seminal paper by Einstein, Podolsky 
and Rosen. Nowadays, we recognize the importance of entanglement as a valuable 
resource for quantum information processing and communication. Its usefulness has 
since been demonstrated clearly in processes like quantum teleportationOE], quantum 
computation^, and quantum cryptography [5]. 

However, concepts such as entanglement and its implications concerning the 
non-existence of a local realism in quantum mechanics have a more fundamental 
role in quantum mechanics. The issue of "locality" as well as notion of quantum 
measurements has given rise to some of the recent and modern interpretations of 
quantum mechanics as well as a better understanding of quantum phenomena H;. It is 
also amidst all these theoretical constructs that Bell proposed an inequality that could 
rule out the hidden variable description of quantum mechanics 7 . Since then, several 
variants of Bell inequalities, some of which were more amenable for experimental 
investigations, have been derived for two-body correlation functions to investigate the 
existence of local realism^. 

Recently there has been much work on the implementation of quantum processing 
on solid state devices. In this paper, we study the thermal states in a system of 
interaction spins and investigate its quantum "nonlocality" . An interesting type of 
entanglement, thermal entanglement, was studied in the context of the Heisenberg 
XXXpiTH]. XXQ31, and XXZQ21 models. The Heisenberg model has been shown to 
have a potential candidate as a model for spin-spin interaction in a solid state quantum 
computer|13|. Being the large Coulomb repulsion limit of the Hubbard model, it has 
been partially realized in quantum dots^2|, nuclear spins^Jj ana - optical lattices[T5]. 
In a recent work, Imamoglu et al have realized quantum information processing 
using quantum dot spins and cavity QED, and obtained an effective interaction 
Hamiltonian based on the XY spin chain between two quantum dots. The effective 
Hamiltonian was shown to be capable of constructing the Controlled-Not gate|16j. 
The XY Hamiltonian is given by 



where 5* = a % /2(i = x, y, z) and a 1 are Pauli operators. When J\ = J2, the XY model 
becomes XX model. In the XY model, the interaction strength between neighboring 
sites is usually assumed to be independent of the sites. In most solid state models 
however, the inter-site coupling strength is site dependent. In this paper we consider 
an extended quantum XX model in which the interaction strength assumes a particular 
site dependent form. 

This paper is organized as follows. In Sec[2 solutions of the extended XX model 
for 4 particles are given. In Sec|3 we construct the temperature dependent correlation 
functions in terms of thermal equilibrium state and investigate the violation of Bell 
inequality for the thermal state. The threshold temperature is given. We also point 
out that the eigenstates of the extended XX model do not realize maximal violation 
of Bell inequality. Effect of external magnetic field is discussed in Sec|3] and we end 
with some discussions in the final section. 
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2. Solution of the extended XX model 

The extended XX Heisenberg model is described by the Hamiltonian 

N-l 

n=l 
N-l 

= ^2 ^n,n+l(o-+CT~ +1 +<7~<7+ +1 ) (2) 
n=l 

where J n ,n+i — \J n(N — n) is the coupling strength between lattices n and n + 1. 
Obviously, the Hamiltonian H describes a nearest-neighbor interaction spin chain. 
Interestingly such a Hamiltonian has been shown to be useful for perfect state transfer 
in quantum spin networks |17j . The Hamiltonian H possesses 2 N complete and 
orthonormal eigenstates. 

When spin chains are subjected to environmental disturbance, they inevitably 
become thermal equilibrium states. The state of a system at finite temperature 
T is given by the Gibb's density operator p(T) = cxp(— H/kT)/Z, where Z = 
Tr[exp(— H/kT)] is the partition function, H is the system Hamiltonian and k is 
the Boltzmann constant, which is set to unity for convenience in this paper. At high 
temperature, the thermal state becomes maximally mixed and do not violate Bell 
inequalities of any kind. It is therefore interesting to consider the critical temperature 
at which a Bell inequality will be violated. For a two-qubit system, we have the 
original Bell inequality. For arbitrary number of qubits, we have the Zukowski-Brukner 
inequality [B]- 

Unfortunately it is not possible to test Zukowski-Brukner inequality for three 
qubits in this case since the correlation functions denned below are zero. Therefore, 
in this paper, we first focus on the next non-trivial case of a 4-qubit system and test 
the violation of local realistic description using the Zukowski-Brukner inequality. The 
extension to arbitrary number of sites, albeit complicating, can also be done in the 
same manner. The Hamiltonian has sixteen eigenvalues 

Eq = E 7 = E s = Ei 5 = 0, 

Ez = E13 = — 1, E± = En = 1, 

Eq = —2, Eq = 2, 

E\ = En = —3, E 2 = E12 = 3, 

E 5 = -4, £ 10 =4. (3) 

The corresponding eigenstates {|</>o), \4>i)i ' ' ' l^is)} can De computed easily and can 
be found in appendix |Appcndix A| The above eigenvalues and eigenstates completely 
determine the thermal states. The density operator p(T) at the temperature T can 
be written as 

fj.=Q 

where = 1/T and the partition function 
Z = Tr(e-P H ) = e-P E » 

= 4 + 4 cosh(3/3) + 4 cosh + 2 cosh(4/3) + 2 cosh(2/3) (5) 
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3. Violation of 4-qubit Bell inequality and the threshold temperature 

To test quantum nonlocality for the state p(T), correlation function Qijki should be 
computed. From the definition of Qijki®, we have 

Qijki = Tr[p(hi ■ (?) ® (hj ■ a) ® (hk ■ a) ® (hi ■ a)] 

15 



1 ^ e ^Tr^) {^(hi ■ a) (hj ■ a) ® • a) ® (n; • a)] 
fj,=a 



ju=0 



where n Q = (sin fto,, 0, cos (9 Q ), a = i,j,k,l. Q^jki is the correlation function for the 
eigenstate \<j>^}, 



Qijki = ^ll^n) (4>n\(ni ■ &) <8> (nj • ct) ® • <?) ® (hi ■ a)} 
For instance, the quantum correlation for the ground state is given by 

%/3 



(7) 



V3 



■ cos ffc cos ft/ sm ft j sin f 



cos( 



: cos ft/ sm fti sm ftfc 



• cos fti cos ft; sm ft,- sm ftfc 



— cos ftj cos ftfc sm fti sm ft/ — 
v/3 



V3 



■ cos t/j cos ft*; sm ftj sm ft/ 



cos ft,/ cos Oj sm ftfc sm ft/ + sm ft,/ sm ftj sm ftfc sm f 



(8) 



Other quantum correlation functions can also be calculated in a similar way. The 
correlation function for the thermal state p(T) are computed using Eq.©. Based on 
the calculated values of Qijki, we construct Bell quantity B 

B = Qim — Q1112 — Q1121 — Q1122 — Q1211 — Q1212 

— Q1221 + Q1222 — Q2111 — Q2112 — Q2121 + Q2122 

— Q2211 + Q2212 + Q2221 + Q2222 (9) 

For a local realistic description, we require —4 < B < 4. In Figure we have 
numerically computed the Bell quantity as a function of temperature. The results 
show that violation of the Bell inequality occurs at T < T — 0.626. We call this 
critical value To the threshold temperature. The maximum value of B for the state 
p(T) approaches 7.917 at temperature close to zero. 

We have also evaluated the Bell quantity B(\<pfj)) in terms of correlation functions 
of each pure state \<p^)- The maximum value of i3(|</> M )) are 

B m ax(\<P»)) =4 for 
6.112 for 



7.917 
5.657 
4.866 
4.060 



for 
for 
for 
for 



^0,15/ 

^1,2,3,4,11,12,13,14) 
05,10) 
06,9) 



(10) 
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Figure 1. For a local realistic description of quantum mechanics, the Bell 
quantity B must necessarily be less than 4. However, the Bell quantity as a 
function of temperature T shows that there is a significant violation of Bell 
inequality at T < 0.626. 



We can explain qualitatively why the maximum value of B for the thermal state 
should be 7.917 by noting that the thermal state p(T) is the linear combination of 
|</v)(</vJ weighted with the factors e~^ E ^. For eigenvalue E 5 = —4, B max (\(f) 5 )) = 
7.917, the power is e 4 ^ and when (3 is large enough, the Bell quantity B is totally 
determined by the contribution of state |0s). Another thing worth noting is that the 
eigenstates of extended XX model do not lead to highest value of B max . We check the 
maximum value of the Bell quantities consist of correlation functions for the following 
three general states 

\<f>') = cos «i| 1000) + sin ai cos a 2 10100) 

+ sinai sina 2 cosa3|0010) + sinai sina 2 sina3|0001) (11) 

\<fi") = cosai|1110) + sinai cosa 2 |1101) 

+ sinai sina2 cosa3|1011) + sinai sina2 sina3|0111) (12) 

\(f>'"} — cosai|1100) + sinai cos a 2 1 1010) + sinai sin a 2 cos as| 1001) 

+ sinai sina 2 sina3 cosa^OHO) + sinai sina 2 sina3 sina4 cosas 
+ sinai sina 2 sin a3 sin a4 sinaslOOll) 

and find that 

BrnaM)) = 6-217 
B m ax{Wo)) = 6.217 
*WK'» = 8-485 

(14) 

for |#>) = 1/2(|1000) + |0100) + |0010) + |0001», K) = 1/2(|1110) + |1101) + |1011) + 
|0111)) and |^") = 1/V6(|1100) + |1010)+|1001)+|0110)+|0101)+|0011)) respectively. 
It is easy to see that the degree of violation of Bell inequality for state \<fi' ) is higher 
than that for the eigenstates \(j>n), (fi = 1, 2, 3, 4) listed in Eq. l|A.2jl . The same results 
also happen for the eigenstates |0 M ), (p, = 11, 12, 13, 14) and \4>^), = 5, 6, 7, 8, 9, 10) 
respectively. We see that among all possible B max , the state \4>'q) yields the largest 
violation. 
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Figure 2. Bell quantity for the cases with magnetic field B 
and 2. 
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4. The effect of external magnetic field 

In this section, we would like to study the effect of magnetic field on the nonlocality 
property of thermal state in a general way, for which the Hamiltonian becomes 

JV-l N 

H' = 2J2 J n,n+i(or+o-- +1 + or-a+ +1 ) + 5 £ a z (15) 

n— 1 n—1 

where 5 is the strength of the magnetic field. It is easy to verify that the eigenstates 
of H' are identical with the ones listed in expression \ A. 2(1 of H, but with different 
eigenvalues. 



E' 



45, — -o-r^-u, i^-oTz-u, 

1 + 25, E' 5 = -4, E' 6 = -2, E' 7 = 0, 
0, E' 9 = 2, E' w = 4, E' u = -3 - 2B, 
J i2 — 3 - 25, E[ 3 = -1 - 2B, E u = 1 - 2B, E[ 5 
and hence, a new correlation function and Bell quantity B' are given 

1 



1 + 25, 



15 

Z' ^ 



ij hi 



-45. (16) 



(17) 



;i2i2 



?2122 



(18) 



& = Qiui — Q1112 — O1121 — Q1122 — Q1211 — ^ 

— Q1221 + O1222 — O2111 — O2112 — O2121 

— O2211 + O2212 + O2221 + O2222 

where Z' = Tr(e~^ H ). Now the violation of Bell inequality depends not only on 
the temperature, but also on external magnetic field. Our numerical calculations are 
shown in Fig [21 

There are five curves corresponding to 5 = 0.1,0.5, 1.0, 1.5, and 2 respectively. 
When 5 = 0.1, the Bell quantity shows a similar variation of the violation of Bell 
inequality as a function of T in the absence of magnetic field. With the increasing 
value of external magnetic field, the maximum value of the Bell quantity approaches 
the value 2 for which the 5 field is about 1.5. The variation of the Bell quantity 
as a function of magnetic field can be explained qualitatively as follows. The p'(T) 
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B 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


To 


0.626 


0.611 


0.556 


0.447 


0.248 


None 


0.122 


0.243 


B 


0.8 


0.9 


1.0 


1.1 


1.2 


1.3 


1.4 


1.5 and above 


To 


0.351 


0.427 


0.467 


0.472 


0.436 


0.343 


0.18 


None 



Table 1. Threshold temperatures for different strengths of the external magnetic 
field. When B = 0.5 and B = 1.5 and above, the values of Bell quantity are no 
greater than 4 at all times. Therefore, no threshold temperatures exist for these 
cases. 

is a different combination of |</yi) (^M compared with p(T). The largest contribution 
of all the states |<^) is determined by the value of B. When B < 0.5, it is the 
eigenstate, \<f>s), which ultimately determines the maximal value of the Bell quantity 
{Bmax = 7.917) since e~ l3E $ = e 4 ^ is the largest power among all the factors. When 
0.5 < B < 1.5, \4>n) takes the place of |0 5 ) with power e (3+2B)/3 and B max = 6.112 
at B = 1.0, for example. When B > 2, e~ l3El5 = e 4B/3 is the one with largest 
contribution and B max = 4. But there are two singular values of B = 0.5 and 1.5. In 
these two cases, B max < 4. The reason for this is that the largest factors of e~^ E ^ are 
e -f3E' 5 = e -PE' lx = e 4/3 for B = o 5j e -pE' 15 = e -f}E' u = e &p for b = 1.5, respectively. 

Thus the Bell quantity is determined principally using a combinations of these two 
elements of Q% H , namely, e^{Q% kl + Q\} kl ) and e^(Q}f kl + Q}j kl ). Note that the 
maximum values of the Bell quantity for the latter two correlation functions are 2.228 
and 2.081 respectively. 

The critical temperatures under different magnetic fields have been found (Table 
H). The variation of To with increasing strengths of B is more complicated. This 
complication arises mainly because the eigenstates contributing to the optimization of 
critical temperatures are different from those needed for the optimization of magnetic 
fields. In the latter case, B max is totally determined by the contribution of state with 
the largest weight or factor for sufficiently large j3. In the former case, depending on the 
value of the external magnetic field, the eigenstates contributing to the optimization 
changes and so the optimization is determined using a combination of the correlation 
functions from different states. In short, the variation of To with B is different from 
that of B max with B. 

5. Conclusion 

In this paper, we consider the extended Heisenberg XX model , modeling the nearest- 
neighbor interaction spin chain. For the 4-qubit extended XX model, it is shown 
that since the correlation functions depend on the temperature and the magnetic 
field, the violation of Bell inequality for the thermal state depends critically on 
these two parameters. The effect of temperature for a local realistic description of 
quantum theory is determined by the threshold value of T below which the thermal 
state violates Bell inequality. The effects of temperature are also studied at different 
strengths of magnetic field. For a fixed temperature, we can find the optimal value 
of the external magnetic field that for the violation of Bell inequalities. Our results 
imply that quantum "nonlocality" could be effectively controlled by magnetic field 
and temperature. We restrict ourselves to the 4-qubit case. However, we could also 
have discussed the violation of Bell inequality for thermal state for 2-qubit and 3-qubit 
cases. For 2-qubit extended XX model, the Bell quantity approaches 2y/2 which is the 



Quantum nonlocality of Heisenberg XX model with Site- dependent Coupling Strengths 



maximal violation of 2-qubit Bell inequality and the corresponding threshold value of 
temperature is T = 0.667 when B = 0. However, for 3-qubit case, the correlation 
function defined by this method is always equal to 0. The violation of Bell inequality 
for arbitrary number of qubit can also be done in the same manner. 
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Appendix A. Eigenstates of the 4-qubit Hamiltonian 



Corresponding to the sixteen eigenvalues of the Hamiltonian 
Eq = Ej = Eg = E\5 — 0, 
E3 — E13 — —1, £4 = £14 = 1, 

Eq = —2, Eg = 2, 

Ei = En = —3, E 2 = E12 = 3, 
E 5 = -4, E w = 4, 

the orthogonal eigenstates are 



(A.l) 



|0000) 

— p(-|1000) + \/3|0100) - V3|0010) + |0001)) 
2v 2 

— p(|1000) + \/3|0100) + V3|0010) + |0001)) 
2v 2 

^L(|1000) - 4=|0100) - 4=10010) + 10001)) 

^L(-|1000) - 4=10100) + 4=|0010) + 10001)) 

-(|1100) - 2|1010) + a/311001) + \/3j0110) - 2|0101) + |0011>) 
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06> 
08> 



in) 

f>12) 

iia> 

iu> 
ii5> 



-(— 11100) + |1010) - |0101) + |0011)) 
^(11100) --||1001) + |0011)) 



_A=(-^|iioo) - ||iooi> + |oiio) - ^|oon» 



0s>) - 2 (_|U00> ~ |1010> + |0101> + |0011>) 

1O ) = 7G 1100) + 2|1010) + V3|1001) + V3|0110) + 2|0101) + |0011» 



^=(-|1110) + V3|1101) - V3|1011) + |0111)) 
^=(|1110) + V3|1101) + V3|1011) + |0111)) 

^=(|iiio) - ^liioi) - ^lioii) + loin)) 

^=H 1110 > - ^l 1101 > + 71 1 1011 > + 1 0111 )) 
inn) 



(A.2) 



Appendix B. quantum correlation functions for each pure states 



The calculation of the quantum correlation functions is straightforward. In this 
appendix, we list all the correlation functions for each eigenstate of the 4-qubit 
Hamiltonian for easy reference. 
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correlation function 


explicit expression 


^ijkl ~ ^ijkl 


cos 8i cos 6j cos 9k cos S/ 


Qijkl Qijkl 


— cos 6i cos 8j cos 6k cos 9i — cos Sfe cos 9i sin 9i sin 9j 

cos 9j cos 9 1 sin 9i sin — f cos Si cos 9i sin Sj sin 9k 

— j cos 6*j cos 9k sin Si sin 9i + cos Si cos Sfe sin Sj sin 9i 

— cos Si cos 9j sin Sfe sin 9i 


(I 2 - o 12 

^ijkl ~ ^ijkl 


—cos9i cos 6*j cos 9k cos 0; + ^ cos 8k cos S; sin 8i sin Sj 
cos 8j cos 0; sin Si sin 8k + § cos Si cos 0/ sin S., sin 8k 

+j cos flj cos 8k sin 0j sin 8i + cos Si cos Sfe sin S., sin 8i 
cos Si cos 9j sin Sfe sin 0; 


O 3 - D 13 

^ijkl ~ ^ijkl 


— cos 8 t cos S.,- cos 8k cos S/ — cos Sfe cos S; sin 6*j sin 8j 

P\ i 

— -^j 2 cos cos S; sin 9; t sin Sfe + j cos Si cos Sj sin 9j sin Sfe 

+| cos 9j cos Sfe sin 9 - L sin S; — cos Si cos Sfe sin S., sin 9i 

— cos Si cos 9j sin Sfe sin 9i 


Qijkl = Qijkl 


— cos Si cos 8j cos Sfe cos Sj + cos 9k cos 0; sin 9i sin Sj 

— ^ cos 8j cos 0; sin 8i sin 8k~ \ cos Si cos Sj sin Sj sin 8 k 

— | cos 0j cos 8k sin Si sin &i — ^ cos Si cos Sfe sin Sj sin 0; 

cos Si cos 9j sin Sfe sin S; 


Qijkl 


cos 8i cos Sj cos 8k cos Sj + cos 8i cos S/ sin Sj sin 8k 
— cos 0j cos 8k sin Si sin 9i — sin Si sin 9j sin Sfe sin S; 


Qjjki 


cos Si cos Sj cos 8k cos S/ + cos Q. cos s in Q i sin $ fc 
cos Q i cos sin Sj sin S/ + | sin Si sin 9j sin Sfe sin S; 


Qijki 


cos Si cos Sj cos 8k cos Sj + C os S., cos S; sin Si sin Sfe 
cos 9i cos Sfe sin Sj sin S; — | sin Si sin Sj sin Sfe sin S/ 


Qijki 


cos Si cos 9j cos Sfe cos Sj — cos Si cos S/ sin 9j sin Sfe 
+ cos 8j cos Sfc sin Si sin S/ — sin Si sin Sj sin Sfe sin S; 


^ijkl 


cos Si cos 9j cos Sfe cos 9i — cos Sfe cos S; sin Si sin 8j 

— ^ cos Sj cos 8i sin S, sin 8k ~ \ cos Si cos S/ sin 9j sin Sfe 

— i cos Sj cos Sfe sin 9 - L sin S; — cos Si cos Sfe sin Sj sin S; 

— ^ cos Si cos 9j sin Sfe sin S; + sin Si sin 9j sin Sfe sin S; 



